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ABSTRACT  OF  THE  DISSERTATION 


Beam  Steering  of  Electrically  Segmented  Piezo-Ceramic 
Ultrasonic  Transducers  Using  Normal  Mode  Coupling 

by 

Hossein  Eslambolchi 

Doctor  of  Philosophy  In  Electrical  Engineering 
(Applied  Physics) 

University  of  California,  San  Diego,  1984 
Professor  Victor  C.  Anderson,  Chairman 

It  Is  well  known  that  normal  mode  coupling  In  large 
diameter  piezoelectric  plates  causes  serious  difficulties  when 
attempting  to  operate  over  wide  frequency  bands.  As  a 
consequence  transducers  are  commonly  constructed  as  a  mosaic  of 
elemental  resonators,  each  of  which  has  a  predominant  single 
mode  of  mechanical  oscillation  at  the  frequency  of  interest. 
Such  transducer  arrays  may  be  electrically  steered  to  angles 
other  than  normal  by  applying  different  phases  of  driving 


voltages  to  different  elements.  A  continuous  plate  can  also  be 
used  to  steer  a  radiated  beam  using  normal  mode  coupling  in  a 
narrow  band  system.  The  technique  is  to  adjust  the  frequency  of 
the  driving  voltage  to  match  the  travelling  wave  velocity  of  the 
normal  mode  which  possesses  the  desired  spatial  phase 
relationship  across  the  face  of  the  plate.  If  the  electrodes  of 
a  continuous  piezoelectric  plate  are  segmented  so  that  regions 
of  the  plate  can  be  driven  with  different  phases,  adjacent 
normal  modes  corresponding  to  a  sine-cosine  spatial  phasing  can 
be  preferentially  excited  to  generate  a  travelling  wave.  A  true 
travelling  wave  will  suppress  the  strong  mirror  lobe  that  would 
exist  if  a  standing  wave  were  excited.  The  velocity  of  the  wave 
will  be  frequency  dependent  and  the  wave  number  can  be 
controlled  to  generate  a  steered  radiated  beam. 

In  this  research  the  theory  of  normal  mode  steering  will 
be  developed,  the  dispersion  curve  will  be  derived  so  the 
travelling  wave  velocities  can  be  evaluated,  and  the  effect  of 
the  mirror  lobe  due  to  the  reflection  from  the  edge  boundary 
will  be  analyzed.  Finally,  the  theory  will  then  be  verified  by 
measurements  on  an  experimental  normal  mode  transducer  which 
will  be  compared  with  a  companion  staved  or  mechanically 
segmented  transducer. 
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1.  Statement  of  the  Problem 

This  thesis  deals  with  electrically  steering  a 
piezoelectric  plate  transducer  to  generate  an  acoustic  beam  at 
angles  other  than  normal  to  the  plate  while  still  retaining  the 
high  coupling  associated  with  normal  mode  operation. 

The  theoretical  analysis  that  treates  the  dispersion 
relation,  reflections  at  the  boundary,  and  grating  lobe  structure 
will  be  presented. 


2.  Introduction 

In  vibrating  piezoelectric  solids,  many  effects  cause 
small  changes  in  the  natural  frequency  of  the  plate,  corresponding 
to  different  vibrational  modes  with  different  amplitudes.  These 
effects  include  such  things  as  material  stiffness,  water  loading 
of  the  surface,  stiffening  of  the  surface  (due  to  the  presence  of 
a  thin  surface  conducting  sheet),  the  biassing  stress,  strains, 
and  electric  fields  which  may  be  applied. 

When  generating  an  acoustic  beam  at  angles  other  than 
normal  to  the  transducer  plate,  we  are  applying  different  phases 
of  driving  voltage  to  different  sections  of  the  transducer.  This 
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phase  difference  causes  a  small  change  in  the  resonant  frequency 
of  the  plate  which  corresponds  to  a  different  mode  of  mechanical 
oscillation. 

Since  the  plate  is  driven  close  to  its  natural  frequency, 
different  modes  with  different  amplitudes  will  be  excited.  But 
since  we  are  applying  different  phases  across  the  plate  as  we 
steer  the  beam  electrically  from  normal,  we  are  particularly 
looking  for  modes  which  have  the  same  spatial  phase  distribution 
as  the  ones  which  we  are  applying.  Therefore  the  problem  is  to 
find  frequencies  close  to  the  plate's  natural  frequency  which  have 
the  same  spatial  distribution.  Finally,  by  matching  these  with 
the  frequency  of  the  driving  voltage,  the  maximum  response  can  be 
achieved.  The  advantage  of  this  continuous  plate  with  normal  mode 
operation  is  that  the  grating  lobe  would  be  suppressed.  By 
contrast,  if  one  applies  the  phase  distribution  to  a  discrete 
segmented  plate  which  has  no  normal  mode  coupling,  the  grating 
lobe  will  pose  a  major  problem  if  the  segments  are  wider  than  one- 
half  wavelength. 
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3.  Background 

Wave  propagation  in  a  homogeneous  infinite  piezoelectric 
plate  has  been  the  subject  of  extensive  investigation  over  the 
years.  Rayleigh  [1]  gave  the  solution  for  an  isotropic  elastic 
plate  and  Ekstein  [2]  gave  the  solution  for  a  particular 
anisotropic  elastic  plate.  Since  then  several  authors  have 
calculated  the  roots  of  Rayleigh  [3-7]  and  Ekstein  [8-10] 
transcendental  frequency  equations  in  order  to  determine  the 
resulting  dispersion  relations  and  solve  high  and  low  frequency 
vibrational  problems.  These  solutions  and  calculations  are 
applicable  to  purely  elastic  plates  and  to  plates  which  obey  the 
linear  piezoelectric  equations  but  which  have  low  electro¬ 
mechanical  coupling  factors.  However,  certain  polarized 
piezoelectric  plates  such  as  PZT-4,  PZT-5  and  PZT-7  which  obey  the 
linear  piezoelectric  equations  have  high  electromechanical 
coupling  factors;  consequently  the  aforementioned  solutions  are 
inadequate  to  describe  wave  propagation  in  these  plates. 

Tiersten  [11]  solved  the  problem  of  wave  propagation  in  an 
infinite  piezoelectric  plate  belonging  to  the  crystallographic 
class  situated  between  shorted  electrodes  for  two  special 
orientations  of  the  sixfold  axis.  The  solution  was  derived  using 
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the  linear  piezoelectric  equations.  He  showed  that  for  a  given 
frequency  and  wave  number  in  the  propagation  direction  there  are 
three  independent  solutions  which  are  coupled  at  the  traction-free 
boundaries  of  the  plate.  The  dispersion  curve  was  derived  from 
the  resulting  transcendental  equation.  Tiersten's  solution  to 
derivation  of  the  dispersion  curve  was  focused  on  the  unforced 
mechanical  vibration  where  there  is  no  electrical  drive  on  the 
sections  of  the  plate.  This  results  in  no  electrical  stress  in 
the  plate. 

In  later  years  Tiersten  [12]  solved  the  problem  of  the 
wave  propagation  by  assuming  that  the  plate  is  driven  electrically 
with  constant  electric  field  in  the  plate  and  traction-free 
boundaries.  The  dispersion  curve  derived  from  his  analysis  is  not 
too  much  different  from  the  undriven  plate.  His  electrical  drive 
consisted  of  uniform  amplitude  and  no  phase  distribution  in  the 
form  of  a  constant  electric  field.  He  investigated  different 
modes  of  wave  propagation  to  determine  which  modes  can  be  excited 
under  electrical  stress.  Both  symmetric  and  antisymmetric  waves 
were  analyzed.  The  different  propagational  modes  contribute 
different  spatial  phase  distributions  in  the  plate  with  different 
acoustical  amplitude  distributions. 
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Bleustein  [13]  approached  the  problem  of  wave  propagation 
in  an  infinite  plate  by  assuming  that  there  are  perfectly  thin 
electrodes  coated  on  the  surface  of  the  plate  and  that  the 
electrodes  are  short  circuited.  His  solution  results  in  zero 
potential  on  the  two  surfaces  of  the  plate  —  a  boundary  condition 
which  is  different  from  that  which  has  been  investigated  in  this 
thesis. 

This  thesis  evaluates  the  dispersion  curve  for  an  infinite 
piezoelectric  plate  by  assuming  pressure  release  boundaries  and  a 
constant  voltage  but  different  phases  to  different  sections  of  the 
plate.  The  theory  will  be  analyzed  in  Chapter  II  and  verified 
against  the  experimental  data  in  Chapter  III. 


4.  Approach 

Much  research  has  been  devoted  to  generating  an  acoustic 
beam  at  angles  other  than  normal  to  the  transducer  plate  using 
piezoelectric  transducers.  Up  until  now,  transducers  have  usually 
been  segmented  by  dicing  the  plate  into  elemental  resonators  which 
have  been  both  electrically  and  acoustically  isolated.  Figure  1 
shows  such  a  segmented  piezoelectric  plate. 


ELECTRICALLY, ACOUSTICALLY  ISOLATED  PLATE 

discrete  TYPE 


For  a  segmented  plate  the  simple  dispersion  equation  can 
be  written  as 


where  Kq  =  wave  number  (discrete  elemental  resonator)  and  Kx,  Ky 
and  K2  are  wave  numbers  along  x,  y,  z,  respectively.  The  equation 
can  be  rewritten  as: 


where  L,  D,  and  T  are  length,  width,  and  thickness  of  an  elemental 
resonator,  n,  m,  and  p  are  integers  (0,1,2,...),  and  Cs  is  the 
velocity  in  the  plate.  Since  the  plate  is  resonant  in  its 
fundamental  mode  of  mechanical  oscillation  in  the  thickness 
direction,  p=l,  m=n=0.  For  such  a  resonator,  there  exists  not 
only  this  fundamental  mode  of  mechanical  oscillation  but  also 
other  higher  order  modes  (n,  m  >  0)  which  are  overtones  in  the 
frequency  response  of  the  plate.  An  array  of  these  resonators  can 
be  electrically  steered  to  angles  other  than  normal  by  driving  the 
individual  elements  at  the  main  natural  frequency,  but  applying 
different  phases  to  the  resonators. 
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An  alternative  approach  that  will  be  demonstrated  here  is 
to  use  a  continuous  plate  to  electrically  steer  a  beam  to  angles 
other  than  normal  by  applying  different  phases  of  driving  voltage 
to  different  electrodes.  If  the  electrodes  of  a  piezoelectric 
plate  are  segmented  such  that  sections  of  the  plate  can  be  driven 
with  different  phases,  then  certain  normal  modes  can  be 
deliberately  excited.  This  can  be  achieved  in  a  continuous  plate 
where  there  is  normal  mode  coupling. 

For  a  continuous  plate  the  dispersion  equation  of  (1)  also 
applies,  only  the  dimension  differs.  The  dimensions  of  a 
continuous  plate  (Figure  2)  are  large  with  respect  to  a  segmented 
discrete  type;  as  a  consequence  there  will  be  normal  modes  which 
are  close  to  the  natural  mode  of  the  plate.  Because  of  the 
coupling  in  the  plate,  there  exist  normal  modes  with  different 
amplitude  and  spatial  phase  distributions  in  the  plate.  In  our 
application  the  plate  is  driven  close  to  its  natural  resonant 
frequency,  and  different  modes  with  different  amplitudes  will  be 
excited.  We  will  particularly  look  for  modes  which  have  the  same 
spatial  phase  distribution  as  that  of  the  impressed  electric 
field.  When  this  match  is  achieved,  we  will  realize  a  maximum 
response  on  the  beam  pattern  of  the  transducers. 


isolated  plate,  continuous  type 
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When  the  frequency  of  the  travelling  wave  which  possesses 
the  desired  spatial  phase  distribution  is  found,  the  reflections 
from  the  edge  boundary  will  create  a  mirror  lobe.  The  creation  of 
the  mirror  lobe  is  due  to  the  standing  wave  that  will  be  formed. 
The  magnitude  of  the  ratio  of  the  travelling  wave  can  be  enhanced 
by  applying  suitable  dissipative  damping  coupled  to  the  edges  of 
the  plate  to  reduce  the  reflections  that  cause  a  standing  wave  to 
be  formed. 

An  active  approach  which  has  been  used  in  this  thesis  is 
to  reduce  the  reflection  of  the  edge  boundary  by  driving  the  last 
few  segments  out  of  phase  and  forcing  a  reduction  in  the  magnitude 
of  the  travelling  wave  by  the  time  it  reaches  the  boundary. 

The  approach  for  a  continuous  plate  also  circumvents  the 
necessity  of  dicing  the  fragile  ceramic  disk,  and  allows  one 
merely  to  use  segmented  electrodes. 


5.  Grating  Lobe  in  Discrete  and  Continuous  Plates 

Discrete  plates  have  disadvantages  due  to  the  grating 


lobe,  which  is  an  aliasing  problem  in  comparison  to  the  continuous 
plate.  In  a  discrete  plate  the  bar  moves  as  a  whole.  The 
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assumption  we  could  make  here  Is  that  the  displacement  across  the 
segment  is  uniform,  so  that  what  we  have  is  a  distribution  that  is 
spatially  quantized,  with  an  accompanying  high  spatial  harmonic 
content  caused  by  the  discontinuities  between  segments.  Because 
of  the  coupling  in  the  plate,  spatial  frequency  content  of  the 
continuous  plate  is  different  in  that  there  are  no 
discontinuities.  This  reduces  the  high  frequency  component  and 
lowers  the  grating  lobe  amplitude  considerably.  The  analysis  for 
the  grating  lobe  will  be  discussed  in  Chapter  II. 


6.  Phase  Quantization 

Phase  quantization  has  been  used  in  this  investigation  in 
order  to  generate  electrical  phase  distribution.  This 
quantization  was  achieved  by  dividing  the  phase  uniformly  between 
zero  and  2*.  Two,  three,  and  four  bit  quantization  was  used  in 
order  to  evaluate  the  effect  of  the  quantization  on  the  beam- 
pattern.  This  quantization  technique  will  be  discussed  in 
Chapter  III. 
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2.0  INTRODUCTION 


The  small  vibrations  of  piezoelectric  bodies  are  governed 
by  the  equation  of  the  linear  theory  of  piezoelectricity.  In 
piezoelectricity  the  quasistatic  electric  field  is  coupled  to  the 
dynamic  mechanical  motion.  To  be  more  specific,  the  equations  of 
linear  elasticity  are  coupled  to  the  charge  equations  of 
electrostatics  by  means  of  piezoelectric  constants. 

Although  there  are  a  reasonable  number  of  papers  on 
piezoelectricity  currently  available  [1-2],  they  have  all  been 
investigations  from  different  practical  points  of  view,  and  use 
the  theory  sporadically,  without  a  systematic  development.  All 
the  existing  papers  discuss  the  piezoelectric  vibrations  of  bodies 
only  in  the  simplest  cases  of  the  thickness  vibrations  and  the 
low-frequency  extensional  and  flexural  vibrations  of  thin  rods. 

In  addition,  considerations  of  large  piezoelectric  coupling  are 
absent  in  papers  on  this  subject  except  for  discussions  of  the 
simplest  case  of  elementary  thickness  vibrations  [3-5]. 

2.1  Dispersion  Curve 

In  this  paper,  the  rigorous  solution  of  the  problem  of 
wave  propagation  in  an  infinite  piezoelectric  plate  is  derived 
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from  the  linear  piezoelectric  equations.  The  surfaces  of  the 
plate  are  pressure  release  and  completely  coated  with  segmented 
electrodes.  No  restriction  on  the  relative  magnitude  of  the 
elastic,  piezoelectric,  or  dielectric  constants  is  Imposed.  The 
solution  will  be  valid  for  materials  with  high  electromechanical 
coupling  factors.  We  will  consider  the  piezoelectric  plate 
polarized  in  the  thickness  direction.  A  transcendental  frequency 
equation  is  derived,  the  roots  of  which  determine  the  modes  of 
propagation  and  the  dispersion  relations  in  an  unbounded  plate. 

Figure  1  shows  a  piezoelectric  plate  which  has  been 
segmented  and  isolated  electrically.  A  rectangular  Cartesian 
coordinate  system  X,,  1*1, ..3  is  chosen  with  X3  =  ±h  defining  the 
faces  of  the  plate  which  is  poled  in  the  thickness  direction.  An 
index  preceded  by  a  comma  denotes  differentiation  with  respect  to 
a  space  coordinate. 
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2.2  General  Equations 

The  system  of  equations  governing  the  behavior  of  this 
plate  consists  of  the  following. 

The  stress  equations  of  motion  [6]: 


The  charge  equation  of  electrostatics: 


The  strain-mechanical  displacement  relations: 


Su  -  1  <UK.«  *  (3) 

The  electric-field  potential  relations: 

EK  =  -*»k  (^) 

The  linear,  piezoelectric  relations  [7-8]: 

Ti  j  "  Ci jK*  SKt  '  eKi  j  EK  (5.1 ) 

Di  =  eiKl  SKt  +  eiKEK  (5.2) 


where,  in  the  above,  Tijf  Uj,  D. ,  S^,  E^  are  the  components  of 
stress,  mechanical  displacement,  electric  displacement,  strain. 
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and  electric  field,  respecti vely;  p  and  ^  are  the  mass  density  and 
the  electric  potential,  respecti  vely ;  eKij*  eiK  are  the 

elastic,  piezoelectric,  and  dielectric  constants,  respectively. 

The  face  of  the  infinite  plate  is  immersed  in  castor  oil 
in  a  surrounding  circular  package  and  tested  in  water.  Since  the 
impedence  of  the  water  is  much  smaller  than  the  impedance  of  the 
ceramic,  we  can  ignore  the  medium  loading  and  basically  assume 
that  we  have  a  pressure  release  boundary.  This  assumption  is  not 
as  valid  at  the  resonant  frequency  of  an  entitled  normal  mode. 

The  faces  of  the  infinite  plate  are  coated  with  segmented 
electrodes  so  it  can  be  steered  electrically.  Since  the 
electrodes  are  assumed  to  be  infinitesimally  thin,  we  can  ignore 
all  possible  mechanical  effects. 

Tjj  =  0  at  =  ±h  (6) 

i  (Kx,  Si  n  8  -  u>t) 

♦  *  RE  [e  ]  at  x3  =  +h  (7) 

where  0Q  is  the  steering  angle,  w  is  the  frequency  of  the  driving 
distribution,  t  is  the  time,  and  ^  is  the  driving  potential. 

For  the  plate  which  has  been  used  in  this  thesis,  lead- 
zirconate,  the  piezoelectric  stress  matrix  can  be  written  as 
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0  0  0  0 
0  0  0  0 

e31  e31  e33  ° 


0  0 
0  0 
0  0 


(8a) 


&2\  states  that  the  plate  is  driven  along  its  thickness  direction 
and  the  resulting  strain  is  along  xj  direction.  Since  the 
material  property  is  the  same  along  the  x1,x?  direction  but  is 
different  along  the  X3  direction,  633  would  be  different  from  both 
e31  and  e32* 

The  dielectric  matrix  can  be  written  as 


e .  . 
ij 


en  0 

0 

0 


'll 


(8b) 


where  again,  because  of  material  symmetry  along 

xl»x2»  en  and  e22  are  the  same  and  e33  dl^^erent. 
The  elastic  constants  are: 


where  p  or  q  =  1,2, 3, 4, 5, 6  may  be  replaced  by  11,22,33,23  or 
32, 31, or  13,12  or  21,  respectively. 

Now  within  the  plate  we  could  write  the  stress  equation  of 
motion  as  follows: 


PUj  = 


pu. 


pu- 


«Tn  6T, 


6T 
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6x? 

Tx.+  “ 

Tx, 

1 

2 

3 

6T21  . 

6T22  . 

6T23 

‘*i 

5x2 

6X3 

6T„, 

6T„„ 

6T„„ 

31  . 

32 

33 

6X1 

"^2 

(9: 


and  by  Maxwell's  electrostatic  relation: 


i-rvvwvv  i"». 


-V  -V-  V'  «V 


Now  by  substitution  of  (3),  (4)  into  (13),  Eq.  (14)  results: 


T11  =  C11  Ul,l  +  C12  U2,2  +  C13  U3,3  +  e31^,3 


T22  =  C12  Ul,l  +  C11  U2,2  +  C13  U3,3  +  e31  *,3 


T33  =  C13  Ul,l  +  C13  U2,2  +  C33  U3,3  +  e33'p,3 


T23  =  C44  <U3.2  +  U2,3> 


T13  -  C44  <U3,1  +  Ul,3> 


T66  =  C66^U2,1  +  Ul,2> 


and  Eq.  (12)  becomes 


Di  =  -  £n  *>i 


d2  ’  ‘  'll*. 2 


°3  “  e31  Ul,l  *  e31  U2,2  *  e33  U3,3  ‘  c33  *,3 


,*  .’*y*  >  /•  .**  .**  .*»  «.*■  i'*  »'*  .*•  .**  v*  ■  »■  *  .*  h'V*  -** »’» 
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Substituting  Eqs.  (14),  (15)  Into  the  equations  of  motion  and 
electrostatic  relations  of  (2),  Eq.  (17)  is  achieved: 


C11  Ul,ll  +  C12  U2,21  +  C13  U3,31  +  e31  **31  +  C66  U2,12 


+  C66  Ul,22  +  C44  U3,13  +  C44  Ul,33  =  ^1 


C66  U2,ll  +  C66  Ul,21  +  C12  Ul,12  +  C11  U2,22  +  C13  U3,32  +  e31  *,32 


+  C44  U3,23  +  C44  U2,33  =  pU2 


C44  U3,ll  +  C44  Ul,31  +  C44  U3,22  +  C44  U2,32  +  C13  Ul,13 


*13  U2,23  +  C33  U3,33  +  e33  *,33  *  PU3 
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"ell  ^,11  ’  C11  *,22  +  e31  Ul,13  +  e31  U2,23  +  e33  U3,33  "  E33  *,33  =  0 

(17) 

and  in  our  new  notation,  the  boundary  condition  becomes: 

T33  =  T32  =  T31  =  0  at  x3  =  ±h  (18) 

Since  the  plate  is  driven  along  its  length  (x^  direction) 
we  could  assume  that  the  mechanical  displacement  is  constant  along 
the  X£  direction. 

Since  the  plate  is  forced  electrically  into  vibration,  the 
total  potential  in  the  plate  can  be  written  as: 

♦  (xlt  x3,t)  =  [’J'h(x1,x3)  +  *ps(x1»x3)]  e1u>t  (19a) 

and  =  0  at  x3  =  +h  (19b) 

where  is  the  homogeneous  steady  state  solution  and  \l>  is  the 

n  p5 

particular  solution  due  to  the  applied  potential  between  the  two 
surfaces  of  the  plate. 


Now  we  need  to  solve  the  homogeneous  solution  where 
(xl»x3)  =  0  at  x3  -h* 

PU1  =  C11  ul,ll  +  c13  u3,13  +  e31*H,31  +c44^u3,13  +  ul,33^ 

pU3  =  C44U3, 11  +  C44U1,31  +  C13U1,13  +  C33U3,33  +  e33*H,33  (20> 

‘  G114'H,  1 1  +  e31ul,  13  +  e33u 3, 33  “  e33*H#33  =  0 
Let's  consider  a  solution  of  (20)  as: 

Uj  =  K1  cosPx3  sin(Lx1  -  u>t) 

u3  =  K3  sinPx3cos(Lx1  -  mt)  (21) 

i|>H  *  B  sinPx3  cos(Lx1  -  wt) 

where  K}, K3.fi  are  constants  and  P.L  are  wave  number  along 
thickness  and  propagation  direction,  respectively.  Notice  the 
form  of  Eq.  (21)  which  is  a  standing  wave  along  the  thickness 
direction  and  a  travelling  wave  along  the  propagation  direction. 
Eq.  (21)  satisfies  (20)  if  the  following  holds: 


w*.*. 
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(^i^L  +  C44P  -  pui  )  Kj  +  (Cjj  +  C44)  LPKj  +  e^^LpB  =  0 

( c44  +  C13)  LPKj  +  (c44L2  +  C33P2  -  pco2)K3  +  e33P2B  =  0  (22) 

e31LPK l  +  e33p2,(3  ‘  (£11l2  +  E33p2^B  =  0 

This  system  of  linear,  homogeneous  equations  in  Kj,  k3,  B 
yields  nontrivial  solutions  when  the  determinant  of  the 
coefficients  of  Kj.Kj.b  vanishes;  i.e.,  when 

( c  1 1 L  +c44^>  “P10  )  ^C13+C44^P  e31^"P 

(C44+C13)LP  (c44L2+c33p2-pa;2)  e33P2  =  0  (23) 

e31  LP  e33P  -^e11L  +e33P  ^ 

Eq.  (23)  is  quadratic  in  u>2,  but  cubic  in  L2,P2.  Therefore,  for  a 
given  w  and  a  propagation  wavenumber  L,  there  are  three  thickness 
wavenumbers  (P^^,  P^2),  p(^)),  each  of  which  yields  an 
independent  solution  and  amplitude  ratios  when  substituted  in 
(22).  The  amplitude  ratios  will  be  designated  by 
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(K 


^  :  B^)  =  :  B2^) 


(24a) 


where 


I 


£ 

.*4 

C- 

i 


6j(i)  •  [(cul2«44P(i)2-p»2)  *  (c44L2*c33P(1)2-pU2)]  -  [e33P(f)2*e33P(i)2] 


B,(1)  «  [e31LP(,)*(c44l2«„P(i)Z-P»2)]  -  l(c44«,3)LP(1)*e33P(1)2] 


,,2 


83(i)  =  [e33P(i)2*e31LP(1)]  *  [(c44«13)LP(i)*(EuL2*e33P(1)  )] 


Eq.  (23)  clearly  shows  that  the  piezoelectric  constants 
e^p  couple  the  quasistatic  electric  solution  of  (17)  to  the 
dynamic  mechanical  motion  of  (20).  These  two  solutions  uncouple 
only  if  all  the  e-jp  vanish.  This  coupling  of  a  quasistatic 
phenomenon  to  a  dynamic  phenomenon  states  that  a  given  frequency 
and  propagation  wavenumber  result  in  three  thickness  wavenumbers. 

One  solution  of  the  differential  equation  is  insufficient 
to  satisfy  the  boundary  conditions;  all  three  solutions  are 
required.  Hence,  we  must  take: 


* 


-T  O  ^ i 


•V A 


>  "V 
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u1  »  sinCLxj  -  iot)  ^  A^^^cosP^x-j 


»  cos(Lxj  -  iot)  J  sinP^x^ 


(25) 


<j<  *  cos(Lx,  -  ait)  Y  A^&«^  sinP^x.. 

1  1-1  *■  * 


as  solutions  of  the  problem. 

Now  we  need  also  to  consider  the  boundary  condition  in 
order  to  solve  the  problem  completely.  Substituting  from  (25) 
into  the  boundary  conditions  (18),  (19b),  and  using  the  required 
relations  in  (14),  (15),  one  obtains: 


T33  *  C13ul,l  +  C33U3,3  +  e33'*',3  =  0 


T13  =  C44*u3,l  +  Ul,3*  =  0  at  X3  =  th 


(26) 


T23  =  C44(u3,2  *  U2,3)  =  0 


«PH  =  o 


and  from  Eq.  (26),  Eq.  (27)  is  achieved: 
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■" 

+ 


I 


I 


T,(i)  sinP(i)h  =  0 


T  A(i)  T?(i)  cosP(i)h  =  0 

frl  2 


(27) 


3 

x 

i  =1 


Y  ^  sinP^  ^h  =  0 


where 


Tl(1)  =  c44  ^3(1)l  +  01(1)  p0)) 


(28) 


T2{1)  C13  ei(1)  L  +  c33  e30)  P(1)  +e33e2(1)  ^ 


(i) 


Eq.  (27)  consists  of  a  system  of  linear,  homogeneous  algebraic 
equations  in  the  A^).  This  system  yields  a  non-trivial  solution 
when  the  determinant  of  the  coefficients  of  the  A^ )  vanishes; 
i  .e.,  when 


D(L,u)  = 


(1)sinP(1^h 

T^^sinP^h 

T^^sinP^ 

^cosP^h 

T2(2)cosP(2)h 

T2(3)cosP(3)h 

^slnP^h 

e,^sinP^h 

6  (3)sinP(3)h 

=  0 


(29) 
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Eq.  (29)  is  a  transcendental  equation,  the  roots  of  which 
enable  the  determination  of  the  dispersion  relation  for  this 
piezoelectric  plate.  Eq.  (29)  contains  an  infinite  number  of 
roots  hn,  each  of  which  determines  a  point  on  the  dispersion 
spectrum  and  yields  amplitude  ratios  (A^1)  :  A^)  :  n(3))  when 
substituted  in  (27). 

Now  for  the  purpose  of  calculating  the  dispersion  curve, 
it  is  convenient  to  write  the  pertinent  equations  involved  in  the 
calculation  in  terms  of  certain  dimensionless  quantities.  The 
pertinent  equations  are  (22),  (23),  (27),  (28),  and  (29).  Let's 
assume: 


Now  the  pertinent  equations  become 


T1  >-J44<Vi)"  +  Sl“(1)> 

4.-33  V1’  “0> 


and  finally,  from  the  determinant  of  (33) 


X(1)  cot  j  a(1)  +  X(2)  COt  ~  a(2)  +  X(3)  cot  j  a(3)  =  0 


where 


x'1’  -  V1’  (V3)  e?(2i  -tV2)  8V3>) 

x(2)  „  -[^(2)  (Tjl1)  8?(3)  .  Tj^3'  B^1') 

X<3>  .  ?  <3>  C^(2>  ;2<»  -  V!)  82(2)) 


The  relations  between  n  and  n,  (i.e.,  the  dispersion 
relation),  can  now  be  evaluated.  Several  methods  of  obtaining  the 
roots  of  the  transcendental  equation  have  been  suggested  by 
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several  authors  [9-10].  But  a  very  simple,  straightforward 
procedure  for  calculating  the  dispersion  relation  is  by  selecting 
a  value  of  n.  The  further  choice  of  ft  permits  the  calculation  of 
the  from  (32)  and  the  §  ^  from  (31)  so  that  the 

and  may  be  computed  from  (34).  The  may  now  be 

determined  from  (36).  If  the  values  thus  determined  satisfy 
Eq.  (35),  the  selected  values  of  n  and  ft  constitute  a  point  on  the 
dispersion  curve.  If  Eq.  (35)  is  not  satisfied,  we  repeat  the 
calculation  for  different  values  of  n  and  ft  until  Eq.  (35)  is 
satisfied  and  a  root  has  been  obtained.  When  a  sufficient  number 
of  values  of  n  and  ft  satisfying  Eq.  (35)  has  been  obtained,  the 
dispersion  curve  can  be  plotted. 

The  starting  point  on  the  dispersion  point  is  critical  in 
determination  of  the  dispersion  curve.  The  solution  at  an 
infinite  wavelength  can  be  obtained  simply  by  setting  L  =  0. 
However,  the  procedure  is  not  quite  as  straightforward.  When  L  is 
set  to  0,  Eqs.  (22)  and  (23)  becomes,  respectively, 

(c44P2  -  pw2)  Kj  =0 

(c33P2  -  pa;2)  K3  +  e33P2B  *  0  (37) 

e33P  K3  "  e33P  B  =  0 

and 


(C44P-pO 


2  2 
C33P  -pu) 


e33P 


"  e33P 


=  0  (38) 


From  Eq.  (38),  the  three  are  given  by 


(P(1),  P(2),P(3))  =  [«  (P/c44)1/2,  0,  u(p/c33)1/2]  (39) 

-  -  2 

where  c33  =  c33  +  e33  /e33  •  The  substitution  of  the  three  P^  ) 
successively  into  the  algebraic  equation  of  (37)  yields 


(6 


(i)  B  (i)  -e  (i)1 

Eq.  (27)  now  becomes: 


A 


(1) 


sinP 


0 


*,2V2) 


cosP^h  +  (c33+e33/£33  cosP 

A^  sinP^h  +  A^  33^ G33 ^  sinP 


<3>h 

<3>h 


0  (40) 


0 


The  last  two  equations  (39,40)  show  that  as  P(2)  ♦  o, 
c{2)  »  since  sinP^h  and  cosP^^h  cannot  equal  zero 

simultaneously.  However,  a(2)  and  p(2)  always  occur  as  an 
indeterminate  product.  Hence,  Eq.  (40)  may  be  written  as: 


A(1)  c44P(1)  sinP^h  =  0 


K  e33  +  A(3)  c33P(3^  cosP^3^h  =  0 


+  A(3)(e33/e33)  sinP^h  =  0 


(21  (21 

where  K  =  Av  'Pv  Eq.  (40)  yields  a  nontrivial  solution  when 


c44P(1)sinP(1)h 


holds. 


e33  c33P^3^cosP^3^h 

h  (e33/e33)s^nP^h 


Eq.  (42)  yields  two  transcendental  equations,  each  of  which 
determines  a  set  of  roots.  The  two  equations  may  be  written  as 
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sinP(1)h  =  0  ,  tanP(3)h  =  P(3)h/Kt2  (43) 

where  *t  =  833) 

Thus  we  see  that  the  first  transcendental  equation  yields 
thickness  frequencies  and  modes  which  are  identical  with  the 
purely  elastic  modes.  We  also  see  that  the  second  transcendental 
equation  yields  thickness  frequencies  and  modes  with  piezoelectric 
coupling.  Thus  for  this  piezoelectric  plate  the  wavelengths  of 
overtone  resonances  are  not  integral  fractions  of  the  fundamental; 
as  a  consequence  the  resonant  frequencies  of  overtone  modes  are 
not  integral  multiples  of  the  fundamental.  The  deviation  from  the 
integral  multiple  relationship  depends  on  the  electromechanical 
coupling  factor  only. 

The  analytical  approach  that  has  just  been  described  here 
was  done  on  the  piezoelectric  plate  (lead-zi rconate)  and  the 
dispersion  curve  was  evaluated.  Figure  2  shows  the  dispersion 
curve  for  such  a  plate.  The  abscissa  is  the  normalized 
(dimensionless)  wavenumber  and  the  ordinate  is  the  normalized 
(dimensionless)  frequency. 


Theoretical  dispersion  curve  (normalized  wavenumber  vs.  normalized  frequency) 


2.3  Electrical  Steerln 


The  evaluation  of  the  dispersion  curve  (frequency  vs. 
wavenumber)  can  now  be  matched  to  the  wavenumber  of  the  driving 
distribution  so  the  beam  can  be  steered  electrically. 

By  Snell's  Law: 

Kq  sineo  =  Kp(«)  (44) 

where  K  =  ,  C  is  the  speed  of  sound  in  the  water,  w  is  the 

0  C  *  a)  r 

b) 

angular  frequency,  and  0q  is  the  steering  angle.  By  solving  (44) 
the  corresponding  frequency  vs.  the  steering  angle  for  such  a 
plate  c,.n  now  be  calculated  from 

K  sine  -  K  (w)  =  0  (45) 

o  o  P 

Figure  3  shows  the  theoretical  curve  for  such  a  plate 
which  has  been  polarized  along  Its  thickness  direction.  The  shape 
of  the  curve  is  similar  to  the  dispersion  curve  of  Figure  2.  The 
asbscissa  is  the  steering  angle  in  degrees  and  the  ordinate,  the 
frequency  in  KHZ. 
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Fig.  3.  Theoretical  dispersion  curve  (steering  angle  vs.  frequency). 


For  an  electrically  segmented  and  acoustically  isolated 
plate  the  acoustic  amplitude  distribution  of  the  plate  is  uniform 
since  the  segments  move  as  a  whole,  the  assumption  we  could  make 
here  is  that  what  we  have  is  a  spatially  quantized  amplitude 
distribution,  with  accompanying  spatial  harmonic  content  caused  by 
the  discontinuities  between  segments. 

Figure  4  shows  a  discrete  linear,  equally  spaced  array  in 
which  the  beam  pattern  function  can  be  derived  as: 

V  =  £  R  cos(u)t  +  mu)  (46) 

m 

where  Rm  is  the  individual  voltage  of  each  array: 


=  tRo 


+  R, 


iu 


i . .+  R 


imu 


me 


+  R 


n-1 


e(n-l)iuj 


where  u  =  KD  sine 

e  *  physical  angle 
D  =  spacing  between  the  elements 


(47) 


K  =  wavenumber 


Fig.  4.  Theoretical  spatial  filtering  function. 


Multiplying  (47)  by  e1u  and  subtracting  from  (47)  results  in: 


V 


a1nU  1 

(1 — zA) 

K  1U  .  ' 

e  -  1 


sin(nu/2) 

sinTu/2) 


(48a) 


Now  the  beam  pattern  V  must  be  multiplied  by  the  directivity  of 
each  individual  element.  That  is, 


V 


sin(nU/2 

sin(U/2 


(48b) 


where  the  second  term  in  (48)  is  the  directivity  of  each  elements. 
Here  if  the  beam  is  steered  off  normal  [11-14] 


u  =  KD(sin6  -  sineQ)  (48c) 

For  the  discrete  segmented  plate,  if  the  path  difference  between 
the  segments  reaches  2tt ,  the  pattern  repeats  itself;  that  is, 

KD(sin6  -  sin90)  =  2mt  n  =  (0,1,2,...)  (48d) 

so  there  exists  some  spatial  harmonic  content  of  the  resonant 
wavenumber  that  gives  rise  to  the  grating  lobe  in  a  discrete 


For  the  continuous  plate,  the  segments'  motion  Is 
different  than  the  one  for  an  acoustically  isolated  plate,  since 
there  are  no  discontinuities  in  the  acoustic  amplitude 
distribution  because  of  the  coupling  in  the  continuous  plate. 
This  amplitude  distribution  has  a  lower  spatial  harmonic  content 
which  reduces  the  grating  lobe  considerably  when  the  beam  is 
steered  off  normal.  The  plate  in  effect  spatially  filters  the 
high  harmonic  that  is  associated  with  the  electric  field. 

This  spatial  filtering  effect  can  now  be  derived 
theoretically  as  a  function  of  a  normalized  wavenumber  in  the 
plate. 

From  (19) 

We  have  already  solved  the  homogeneous  solution  of  the 
problem  when  the  dispersion  curve  was  derived.  Now  we  need  to 
find  the  particular  solution  due  to  the  applied  electric  field. 
From  the  particular  solution  we  obtain  the  acoustic  amplitude 
distribution  as  a  function  of  the  wavenumber  in  the  plate  for  a 
fixed  frequency. 
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Since  the  plate  is  driven  externally  with  potential  of 
Eq.  (7)  at  the  faces  of  the  plate,  there  will  be  a  corresponding 
stress  and  electric  field  due  to  the  applied  drive.  Eq.  (5.1)  and 
(5.2)  can  now  be  written  as: 


T.  .  , 

ij.i 


D.  .  +  0.  .  =  0 

i,i  i.i 


(49a) 


where  T-j j ,  uj,  are  the  components  of  stress,  mechanical 
displacement,  and  electric  displacement,  respectively.  The 
quantities  D*T  represent  the  electrical  term  due  to  the  applied 
potential.  The  associated  boundary  conditions  for  a  pressure 
release  boundary  with  zero-normal  component  of  electric 
displacement  is  given  by: 


N. (T.  .)  =  0 

i  ij 


Ni(D.+0.e)  =  0 


(49b) 


The  previously  derived  homogeneous  solution  yielded  the 
functions  (u^1,  u™,  'J'™),  which  satisfy  the  eqution  of  motion  (49a) 
and  boundary  conditions  (49b).  Now,  if  we  assume  that  the 
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solution  due  to  the  applied  field  is  (u^,  Uj,  4>)»  these  functions 
must  also  satisfy  (49a). 

Now  let's  multiply  (49a)  by  U1*1  once  and  then  by  Ui  and 

J  v 

subtract  and  integrate  over  the  volume: 


l  uj  -  +  ““'“j*  ^ dV  =  0 


Now  by  the  divergence  theorem: 


<“»  -  “Z>  /  0  ft  dV  -  I  V0-*"  -  ft  ds 


Since  the  homogeneous  solution  yielded  zero-normal  component  of 
stress,  mechanical  displacement  and  electrical  displacement  at  the 
boundary,  Eq.  (51)  can  be  rewritten  as: 


<“m  -  “2>  {  ft  uj  dV 


/  N  (Df*)  dS 
S  1 


Now  we  can  write 


■j-.V.’ 


•b  -  W-  •  N  =  i p  u“  UJ dV  (54 

Eq.  (53)  shows  that  at  w  =  oj  the  displacement  function  goes  to 

-i  >  /  m 

infinity,  and  this  is  usually  the  case  for  the  driven  plate  if 
there  is  no  damping  in  the  plate.  Now  the  displacement  and 
potential  function  due  to  the  applied  field  is  required,  so  first 
we  expand  the  solutions  of  the  problem  in  Fourier  sine  and  cosine 
series: 

U  -  I  V*3>  s1n(Lnxl> 

n=l 

U3  '  %  W  cos<Ln’‘l)  <56> 

n-i 

oo 

t  =  I  B  (x  )  cos  (L  x  ) 
ps  n  3  n  1 

where  U^,  U3  are  the  displacements  along  xj,  X3,  res¬ 
pectively,  »|>  is  the  particular  solution  due  to  the  applied 
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electric  fields,  U  (x_),  V  (x0),  B  (x0)  are  acoustic  amplitude 
n  3  n  3  n  3 

distributions  in  the  plate,  and  Ln  are  the  discrete  eigenvalues  of 
the  wavenumber. 

Substitution  of  (55)  into  the  linear  piezoelectric 
equations  of  (5.1),  (5.2)  gives 

K1  U(X3}  +  (K2  +  K4}  V*3  +  K3B,3  +  K5V*33  =  0  (56a) 

L1V(X3}  •  ^2  +  K4}  U*3  +  L2  V‘33  +  L3B*  33  =  0  <56b> 

M1B(x3)  '  K3  U*3  +  L3v» 33  +  M2B*33  =  0  (56c) 

where 


II 

fH 

-C"llLn 

K2  =  "C"l3Ln 

K3  =  "e31Ln  K4  =  “^44Ln  K5  =  ^44 

4  = 

-C  L2 
H'44Ln 

L2  =  C33 

L3  =  e33 

II 

r-H 

r 

C  L2 
LllLn 

M2  =  "e33 

(57) 

The  sets  of  equations  in  (55)  constitute  the  solution  for 
the  forced  vibration  problem  considered  here. 


Now  multiply  (56b)  by  -K3  and  (56c)  by  (K2+K4)  and  add: 


V, 


33 


x0  x, 

iir) v  +  (1 r ) B* 


'1 


'1 


33 


where 


X1  =  L3  (K2  +  V  -  K3L2 

x2  -  k3  4 

X3  =  K3  4  *  M2(K2  +  K4) 
x4  -  Mx  (K2  ♦  K4) 


(58) 


(59) 


Now  taking  the  derivative  of  the  (56a)  with  respect  to  the  x3 
coordinate  and  multiplying  it  by  (K2+K4),  and  (56b)  by  Kj,  and 
adding: 


U,1  +  (K2  +  V,33  +  K3B,33  +  K5  U* 333 

4V  -  (K2  ♦  k4)  u,3  ♦  4  v,33  +  l3  b,33 


0 

0 
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which  results  in 


V*33Y1  +  VY2  +  B*33  Y3  +  U*333Y4  =  0  <6°) 

where 

Y1  ■  <K2  *  *■/  *  h  L2 

V2  *  L1K1 

Y3  ■  K3  (K2  ♦  *  K,  4 

V4  "  K5  (K2  *  V  (61> 


Also  from  (56c),  let's  take  the  double  derivative  with  respect  to 
X3  coordinate: 


M1  L3  M2 

Kj  B,33  +  Kj  V*3333  +  B*3333 


Substitution  back  into  the  equation  of  (60): 


(6: 


V,3333^  +  V*  33Y1  +  VY2  +  B,33^Y3  + 


M, 


+  B * 3333  ^Y4(K3^  =  0 


Now  let's  take  the  double  derivative  of  (58)  and  substitute  back 
into  (58): 


o  x9x^  x.  x?x.  x?  2 

^ * 3333  =  77  B * 3333  +  (“T"  "  77^  B*33  +  B  +  V  (6 


'1 


X1  1 


ll 


'1 


By  substituting  (64)  into  (63),  Eq.  (65)  results: 


L-x-,  M9  x9x,  x.  x9  _ 

Ciqs7  *  4’ B’3333  *  l^T '  4  *  4 ((K2  * K4>  *  k'L2) 


+  (K3(K2  +  K^)  +  K^L^)]  B»33  +  [ 


X2X4  X4 


*-  A-t 

X1  1 


((K2+  k4)2  +  k1l2)]  b 


X2  2  x2 


+  ^7f)  --^K2  +K4)  +Kl4)  +KlLll  v  =  0 


(6 


in  which  V  can  now  be  defined  as  a  function  of  B: 
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V  ‘  (if)  ® * 3333  *  (if)  b-33  *  (if)  B 


where 

L3x3  M2 

zi  =  ^  +  <7 


X2X3  X4  X3 


IT  '  x7  +  xf  «<K2  +  K4>  +  *l4)  +  (K3<K2  +  K4>  +  K1L3) 

1  1  1 


x2x4  x4 


t—  Ai 

xi  1 


•if((K2*K4>  *K1L2» 


:4  ■  if  (^  *  S'"  *  Kl4)  -  K1L1  -  (if)  ) 


Now  by  substituting  (66)  back  into  (58): 


N1Bi)(VI)(x3)  -  "2Bn(IV!<x3>  *  M3Bn(II,(x3)  *  =  0  (68) 


where 


1  =  ^"^13Ln~>~Ln^31Lnc33^^'^31Ln~e33^13Ln'K~n^^33^13Ln“Ln^31Ln^33^ 
"{e31Ln^  Ln^e31Ln^  "  cllLn  "  ^c13Ln+Ln^  "  4lLn  *33^ 
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+e31Ln  eU(c13Ln+  Ln)  r  ,  2ft  ?  ,  .  *3  .  2r  n 

=~~-=—  — - r— —  -  en  Lnl(-c13f-n“Ln>  '  cULn  c33  J 

e33^c13Ln+Ln)  e31Lnc33 


^31Ln+£33^_^13Ln‘Ln^^31Ln^  "  ^ULn“  L  .  ,  V  -  ,  -  1 

(-c,0L  -L  )  +  e0iL  c,0 
'  13  n  r  31  n  33 


®31  Ln  (l  -  ,  ,  .2-  ,2  r  ^  ,4  ,  e3l  Ln  >2^ 

~  -  — .  "  -'7  ^‘c13Ln'Ln)  “cULn  c33 J“cllLn“  CT-; — T  ^ 

(_c13Ln"Lr\)+e31Lnc33  ( “c13Ln‘Ln)+e31Ln 


,3.  re31Ln  Ell(  c13Ln4‘n)  ellLn  (c13Ln^n  ff  -  .  .  >2  -  .  2- 

(  31  n)  l77TTTJT7^  '  TTTTTTTT^  K^w'-nV  -hiVWJ- 

(~c13Ln  ^n ^  e31Lnc33  (c13Ln^n  )_c31Lnc33 


( _ e31  _ irr/  r  i  j  \2  -  1 2  -  \  -  .4  _ e31  ^n _ 1 

^  ,  J  v  -  ,  -  )i(("c13Ln_Ln)  _cllLn  c33-*"cllLn'  .  -  ,  ,  .  -  J 

( "c13Ln  "Ln  ^  e31Lnc33  *  "c31Ln  ^e31Lnc33 


"  ell  Ln  (c13  Ln  +  Ln) 


The  solution  to  Eq.  (55)  will  be  the  acoustic  amplitude 


distribution  in  the  plate,  the  particular  potential  solution  now 
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must  be  matched  at  the  boundary  of  the  plate  to  complete  the 
solution  of  the  problem.  Eq.  (68)  is  a  sixth-order  ordinary 
differential  equation,  which  has  six  roots,  but  because  of  the 
form  of  the  solution  only  three  roots  need  to  be  considered,  since 
the  other  three  roots  yield  the  same  solution  and  need  not  be 
considered. 

Now  let's  take  the  solution  of  the  form 

■y  ^ 

B  (x3)  =  E  e  3  (70) 

Substitution  in  Eq.  (68)  results  in: 

MjY6  +  M2y4  +  M3y2  +  M4  -  0  (71) 

Since  only  three  roots  need  to  be  considered,  this  results  in: 

MjY3  +  M2y2  +  M3y  ♦  M4  =  0  (72) 

The  roots  of  (72)  must  now  be  introduced  back  into  the  equation  of 
(70): 


B  (x3)  = 


Ei  e 


V3 


i  =1 


(73) 


Also  from  (55) 


*  3  Y.x- 

'PpcUpXo)  =  I  1  E.  e  cos  (Lx.)  (74) 

1  *  n=l  i=l  1  n  1 

which  must  also  satisfy  the  driving  potential  at  the  face  of  the 
plate  (i .e.,  X3  3  +h) 

00 

^(Xj.h)  =  £  B  (h)  cos(L  Xj)  (75) 

n=l  " 

where  Lp  now  must  be  matched  with  the  driving  wavenumber  (i.e., 

K  =  —  sine  ),  where  in  the  above  c  is  the  speed  of  sound  in 
co  w 

OJ 

water,  oj  is  the  driving  frequency,  8q  is  the  steering  angle  and 
Bn(h)  is  the  amplitude  of  the  wave  at  the  surface  evaluated  from 
(73). 

To  obtain  the  spatial  filtering  curve  (i.e.,  acoustic 
amplitude  distribution  versus  normalized  wavenumber),  a  particular 
wavenumber  corresponding  to  specific  steering  angle  was  chosen  and 
the  operating  frequency  was  set  at  the  corresponding  resonant 
frequency.  This  results  in  the  evaluation  of  the  acoustic 
amplitude  distribution  from  (75). 


•  *  *  *  »  *  1  j  '  *  '  ,  *  •  *  »  »  _"*•  , 


59 


The  theoretical  analysis  that  was  just  derived  was 
performed  for  the  parameters  of  the  piezoelectric  material  used  in 
this  thesis  for  a  fixed  normalized  frequency  and  10°  steering. 
Figure  4  shows  such  a  plot.  The  abscissa  is  the  normalized 
wavenumber  and  the  ordinate  is  the  normalized  acoustic  amplitude 
distribution  in  the  plate  truncated  to  unity.  It  is  clear  that  at 
the  normalized  wavenumber  of  10°  steering,  the  acoustic  amplitude 
distribution  goes  to  infinity,  and  as  the  wavenumber  changes,  the 
amplitude  decreases  exponentially.  This  spatial  filtering  curve 
is  the  main  phenomenon  in  the  continuous  piezoelectric  plate  that 
serves  to  reduce  the  grating  lobe. 

Since  the  phases  to  steer  the  beam  in  the  continuous  plate 
have  been  quantized,  the  numerical  analysis  using  fast  Fourier 
transform  was  performed  to  get  the  theoretical  beam  pattern  for  a 
continuous  plate.  The  beam  pattern  for  a  continuous  piezoelectric 
plate  length  of  2L  can  be  written  as: 

V(K-.o)  -  JL  4K'V1ne  (76) 

-L 

where  c(x)  is  the  quantized  phase  necessary  to  steer  the  beam 
electrically.  The  first  term  in  the  integral  is  the  directivity 


pattern,  K'  is  the  wavenumber,  xj  is  the  distance  (i.e.,  -L  <  X}  < 
L)  and  8  is  the  physical  angle.  Taking  the  inverse  Fourier 
transform  of  the  second  term  in  Eq.  (76)  gives; 

D(K)  =  J  lkx  [costc(x)]  -  i  sin[c(x)]]  dx  (77) 

Now  the  spatial  filtering  function  of  Figure  4  must  be  multiplied 
to  obtain  a  new  spatial  distribution: 

Y(K)  =  Z(K)  D(K)  (78) 

where  Z(K)  is  the  spatial  filtering  curve  defined  numerically  from 
(75).  Now  by  taking  the  inverse  Fourier  transform  of  (78),  we 
obtain 


Y(x)  =  /  elKx  Y(K)  dK 

-L 


(79) 


and  the  beam  pattern  can  now  be  evaluated  by  taking  yet  another 
inverse  Fourier  transform  by  using  Eq.  (76): 

V ( K ' ,0)  =  /L  £K'X1  sin9  Y(x)  dx 

-L 


(80) 


Figure  5  shows  the  theoretical  beam  pattern  for  2\  spacing 
where  there  is  no  spatial  weighting  function  and  the  beam  has  been 
steered  to  10°  with  uniform  amplitude  shading.  The  first  grating 
lobe  can  be  seen  to  be  at  about  27°.  Figure  6  shows  the  same  beam 
pattern  identically  except  the  spatial  filtering  function  has  been 
introduced.  The  grating  lobe  can  be  seen  to  be  reduced 
considerably.  In  contrast  to  the  acoustically  isolated  and 
electrically  segmented  plate  the  acoustic  amplitude  distribution 
in  the  plate  is  uniform  since  the  segments  move  as  a  whole,  so 
there  is  no  spatial  filtering  effect  in  a  discrete  segmented 
piezoelectric  plate. 


Theoretical  beam  pattern,  10°  steering,  2X  spacing,  uniform 
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CHAPTER  III 


3.0  Introduction 


In  this  chapter  the  theoretical  analysis  of  Chapter  II 
will  be  verified  against  experimental  data.  The  first  part  of 
this  chapter  provides  a  description  of  the  receiver-transmitter, 
the  physical  condition  in  which  the  transducer  was  tested,  the 
experimental  set-up,  system  design  and  the  phasing  scheme  used  to 
steer  the  beam  electrically  over  a  certain  aperture  range. 

In  the  second  part  of  this  chapter  the  experimental 
measurements  for  the  dispersion  curve,  the  grating  lobe,  and  the 
mirror  lobe  in  the  continuous  plate  and  in  the  discrete  plate  when 
applicable  will  be  analyzed.  The  effect  of  two,  three  and  four 
bit  quantization  on  side-lobe  structure,  the  mirror  lobe  and 
grating  lobe  will  be  compared  against  the  lobes  of  a  segmented 
discrete  type. 

3.1  Recei ver/Transmitter 

In  this  part  the  receiver  and  transmitter  properties,  the 
physical  condition  under  which  the  transducer  was  tested,  the 
instrumentation  system  design  to  obtain  the  data,  and  the  phasing 
scheme  used  to  steer  the  beam  electrically  will  be  discussed. 


3.1.1  Transducer 

Although  the  transducer  used  in  this  thesis  to  provide  the 
experimental  measurement  could  be  used  both  as  a  receiver  and 
transmitter,  it  was  used  as  a  transmitter  only.  The  transducer 
material  used  is  piezoelectric  ceramic  (CHANNEL  INDUSTRIES 
5500).  CHANNEL  5500  has  a  high  sensitivity  and  good  time 
stability.  The  dimensions  of  the  piezoelectric  plate  used  in  this 
experiment  are  2  x  2  x  0.24  inches.  It  is  operated  in  its 
thickness  fundamental  mode  of  mechanical  oscillation  and  is  also 
polarized  along  its  thickness  direction. 

The  numerical  values  of  elastic,  piezoelectric,  mass 
density  and  dielectric  constants  obtained  from  CHANNEL  INDUSTRIES 
for  lead  zirconate  (CHANNELL  5500)  are: 

Cjj  =  12.1x10*°  Newton/meter^  p  =  7.75xl03  kg/m3 

in  ? 

C33  =  11.1x10  Newton/meter  e^/eQ  =  1730  Relative  dielectric 

constant,  free 

in  2 

C^4  =  2.11x10  Newton/meter  e33/e0  =  1700  Relative  dielectric 

constant,  free 

in  2  ? 

C 12  =  7.54x10  Newton/meter  e^  =  15.8  Columb/m 

C,,  =  7.52x10*°  Newton/meter^ 
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For  CHANNEL  type  5500  piezoelectic  material  the  frequency 
constant  of  the  thickness  mode  is  approximately  78  KHZ-INCH. 

Using  this  figure  the  resonant  frequency  of  the  continuous  plate 
was  calculated  to  be  approximately  325  KHZ.  The  actual  maximum 
response  at  zero  degree  steering  was  measured  to  be  at  348  KHZ. 

In  order  to  permit  the  beam  to  be  steered  electrically, 
the  transducer  was  segmented  by  a  diamond  saw  to  isolate  the  thin 
conducting  sheet  electrically.  Twenty  segments  approximately  0.92 
inches  wide  were  cut.  For  the  discrete  type  the  segments  were  cut 
deep  enough  to  isolate  the  segments  both  acoustically  and 
electrically.  Figure  1  shows  the  physical  cross-section  for  such 
a  plate.  The  segments  for  the  continuous  plate  were  isolated 
electrically  by  cutting  only  the  thin  conducting  sheet  with  very 
shallow  depth.  Figure  2  shows  the  cross-section  for  this  plate. 

Both  transducers  were  set  in  polyvinyl  circular 
packages.  The  piezoceramic  elements  were  mounted  only  at  the  top 
and  bottom  edges  of  the  plate,  where  the  acoustical  loading  is  so 
minimal  that  it  can  be  neglected.  Figure  3  shows  the  cross- 
section  of  the  package  which  was  used  in  this  thesis.  Since  the 
plate  segments  will  be  driven  electrically  the  package  was  filled 
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Electrically  acoustically  isolated  plate,  discrete  plate 


Electrically  isolated  plate,  continuous 


with  castor  oil.  The  main  reason  for  choosing  castor  oil  as  an 
insulating  medium  is  that  its  acoustic  impedance  matches  that  of 
the  water. 

Since  the  plate,  which  can  have  a  wide  variation  in 
electrical  impedance,  must  be  driven  electrically  over  long  wire 
length,  there  could  be  additional  phase  shifts  which  could  distort 
the  phasing  scheme  and  could  affect  the  beam  pattern  of  the 
system,  especially  the  side-lobe  structure.  A  solution  to  this 
was  to  shunt  the  segments  with  a  resistor  which,  in  combination 
with  the  1.8  KOHM  reactance  of  the  ceramic  (at  resonance)  yields  a 
characteristic  impedance  of  approximately  100  OHM,  matching  the 
characteristic  impedance  of  the  twisted  pair  wire.  The  cable  used 
is  made  by  the  3M  Company  as  ribbon  cable  type.  It  has  20  twisted 
pairs  which  are  used  for  the  twenty  segments. 

The  beam  width  of  the  transducer  used  in  this  experiment 
was  calculated  to  be  approximately  5.2  degrees  at  the  operating 
frequency  of  325  KHZ  beyond  1  meter  and  5  cm  wide  at  shorter 
ranges. 

The  hydrophone  used  to  obtain  the  acoustical  data  is  an 
encapsulated  1/8x1 /8  inch  ceramic  cylinder  resonant  at  about 
380  KHZ.  It  was  located  approximately  1  meter  away  from  the 
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transducer  so  that  a  far  field  measurement  of  the  transducer  would 
be  valid.  This  distance  was  calculated  to  be  1/2  meter 
from  0  /X,  where  0  is  the  aperture  size  of  the  transducer  and  X  is 
the  wavelength  at  the  operating  frequency  of  325  KHZ. 

The  experiment  was  tested  in  a  2'x3'x4'  polyethelene-lined 
tank.  Figure  4  shows  the  tank  in  which  the  system  was  tested. 
Because  of  the  limited  dimensions  of  the  test  tank,  the  pulse 
width  for  the  signal  was  set  short  enough  to  isolate  the 
reflection  from  the  boundaries. 

3.1.2  Instrumentation  System  Design 

Figure  5  is  a  block  diagram  of  the  instrumentation  system 
used  in  this  thesis  which  consists  of  seven  basic  blocks: 
transducer  driver,  transducer,  hydrophone,  hydrophone  receiver, 
analog  to  digital  converter,  system  motor  drive  and  data 
acquisition  system. 

The  system  was  designed  to  provide  arbitrary  phasing  of 
the  individual  segments.  This  was  set  by  picking  a  minimum 
quantization  interval  of  2n/16.  To  obtain  this  we  multiplied 
the  carrier  frequency  by  16  on  a  frequency  burst  generator 
(WAVETEK)  to  clock  the  sixteen  stage  shift  register  (74LS164).  A 
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Instrumentation  system  design  diagram 
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divide  by  sixteen  counter  output  (74LS163)  was  used  as  the  input 
to  the  shift  register.  This  method  provided  sixteen  phases 
uniformly  between  zero  and  2u.  The  first  eight  phases  out  of  the 
shift  register  were  used  as  inputs  to  a  pair  of  (8x10)  mini -matrix 
selector  switches  and  the  last  eight  to  another  pair  of  (8x10) 
switches.  These  selector  switches  provided  the  selection  of  any 
of  the  sixteen  phases  to  each  of  the  20  individual  segments. 

Twenty  double-pole  single  throw  switches  were  used  to  select 
between  the  first  eight  and  the  last  eight  phases.  A  variable 
gain  amplifier  is  used  to  amplify  the  signal  that  now  must  drive 
the  transducer  segments.  This  variable  gain  amplifier  allows  the 
amplitude  to  be  adjusted  on  each  individual  segment  to  obtain  any 
desired  drive  amplitude  shading  across  the  transducer. 

The  hydrophone  output  was  filtered  to  reduce  the  noise  in 
the  channel  and  increase  the  signal  to  noise  ratio  for 
measurements.  A  full  wave  rectifier  and  a  low-pass  filter  was 
used  to  generate  the  envelope  of  the  acoustic  signal  and  then 
amplified  by  two  stages  of  low  noise  and  high  bandwidth  gain 
(HA4620-8)  to  drive  the  A/D  converter  input. 

Block  five  shows  the  analog  to  digital  converter.  This  is 
a  multiplexed  A/D  converter  controlled  over  the  IEEE-488  bus 
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(Connecticut  Micro  Computer).  The  A/D  converter  has  sixteen 
channels  of  analog  input  multiplexer  which  sequentially  sample  the 
input  data  at  a  110  microsecond  sampling  rate  into  an  8  bit  packed 
data  stream.  It  also  includes  the  time  which  was  used  both  to 
trigger  the  WAVETEK  and  to  reset  the  divide  by  sixteen  counter  and 
shift  register  after  the  analog  data  was  converted  into  digital 
data. 

Block  six  is  the  system  motor  drive  which  is  controlled  by 
the  computer.  This  motor  is  a  SLO-SYN  motor  with  72  steps  per 
revolution.  In  order  to  achieve  higher  resolution  on  the  rotating 
angle,  a  reducing  gear  was  used.  This  motor  was  rotated 
mechanically  by  using  a  two  bit  sequence  that  was  controlled  by 
the  software  in  the  computer  via  its  parallel  printer  port. 

Block  seven  is  the  data  acquisition  system  which  is  the 
main  processor  in  the  experiment.  This  microcomputer  (ACCESS) 
uses  its  IEEE-488  port  to  control  the  WAVETEK  parameters  and  to 
control  and  receive  data  from  the  A/D  converter.  It  also 
generates  the  sequence  needed  to  step  the  motor  via  the  parallel 
printer  port.  The  data  reduction  was  controlled  in  the  software 
and  then  stored  on  a  5  1/4  inch  floppy  disk  (MAXWELL  MD1).  The 
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divide  by  sixteen  counter  output  (74LS163)  was  used  as  the  input 
to  the  shift  register.  This  method  provided  sixteen  phases 
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of  the  sixteen  phases  to  each  of  the  20  individual  segments. 

Twenty  double-pole  single  throw  switches  were  used  to  select 
between  the  first  eight  and  the  last  eight  phases.  A  variable 
gain  amplifier  is  used  to  amplify  the  signal  that  now  must  drive 
the  transducer  segments.  This  variable  gain  amplifier  allows  the 
amplitude  to  be  adjusted  on  each  individual  segment  to  obtain  any 
desired  drive  amplitude  shading  across  the  transducer. 

The  hydrophone  output  was  filtered  to  reduce  the  noise  in 
the  channel  and  increase  the  signal  to  noise  ratio  for 
measurements.  A  full  wave  rectifier  and  a  low-pass  filter  was 
used  to  generate  the  envelope  of  the  acoustic  signal  and  then 
amplified  by  two  stages  of  low  noise  and  high  bandwidth  gain 
(HA4620-8)  to  drive  the  A/D  converter  input. 

Block  five  shows  the  analog  to  digital  converter.  This  is 
a  multiplexed  A/D  converter  controlled  over  the  IEEE -488  bus 


plot  of  the  beam  pattern  was  made  on  the  integral  dot-matrix 
printer  of  the  ACCESS  by  a  software  plotting  package  obtained  from 
Enercomp  Company. 


3.1 .3  Phasing  Scheme 

Quantization  of  the  phase  was  necessary  in  order  to  steer 
the  beam  electrically.  The  required  phases  can  be  calculated  for 
four,  three  and  two  bit  quantization  and  for  different  spacing 
between  the  segments.  The  input  drive  to  the  system  can  be 
written  as: 


V(t)  »  A  COS  ( IDT  -  $) 

where  A  is  the  amplitude  of  the  drive,  id  is  the  driving 
frequency,  t  is  the  time,  and  <|>  is  the  phase.  The  phase  which 
must  now  be  quantized  can  be  written  by  setting  V(t)  equal  to 
zero;  (i.e.,  $  =  idt  =  KC^t)  where  C^  is  the  speed  of  sound  in 
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water.  Since  the  plate  has  been  segmented,  the  total  path  can  be 
written  as: 

x(N)  =  C  t  =  C  NO  sin0 

where  N  is  the  segment  number,  0  is  the  segment  spacing.  This 
results  in: 

4(N)  =  INT[ KN°^-i n— -]  MOD  L 

where  L  is  the  quantization  variable,  which  can  be  sixteen,  eight, 
or  four  depending  on  the  quantization  desired  (4,  3  or  2  bits). 

Figure  6  shows  an  example  of  four,  three  and  two  bit 
quantization  for  10°  steering,  one-half  wavelength  spacing. 

Figure  7  shows  an  example  of  four  bit  quantization  at  10°  steering 
for  different  spacing.  For  IX  and  2x  spacing,  the  segments  were 
electrically  tied  together  to  form  larger  segments. 

3.2  Experimental  Analysis 

In  this  part  of  the  thesis  the  experimental  dispersion 
curve,  grating  lobe  and  mirror  lobe  will  be  discussed  for  both  the 
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Fig.  7  Quantization  of  n/8,  y  X,  X,  2X  spacing 


discrete  and  continuous  types  of  transducers  when  applicable.  The 
effect  of  quantization  on  the  beam  pattern  of  the  two  transducers 
will  be  evaluated. 

3.2.1  Dispersion  Curve 

The  theoretical  analysis  that  was  derived  in  Chapter  II 
was  evaluated  for  the  parameters  of  the  experimental  continuous 
segmented  transducer.  In  order  to  obtain  the  dispersion  curve 
(i.e.,  change  in  velocity  versus  frequency),  a  set  of  phases 
computed  for  a  desired  steering  angle  was  set  up  on  the  segments 
of  the  transducer.  Then  a  sequence  of  beam  patterns  similar  to 
Figures  8  and  9  was  generated  to  identify  the  frequency  of  maximum 
response  out  of  the  transducer.  The  angle  of  the  maximum  response 
was  then  measured  and  compared  against  the  desired  steering 
angle.  If  the  angles  were  matched  then  the  frequency 
corresponding  to  the  maximum  amplitude  constitute  a  point  on  the 
dispersion  curve;  if  the  angle  was  not  matched,  another  frequency 
was  chosen,  another  set  of  phases  computed,  and  another  beam 
pattern  sequence  was  run.  This  Iteration  process  was  repeated  as 
necessary  to  obtain  points  on  a  dispersion  curve. 


Figure  8  and  9  show  the  final  beam  pattern  for  a  5° 
steering  angle  at  one-half  wavelength  segment  spacing.  The 
abscissa  is  the  angle  in  degrees,  the  ordinate  is  the  normalized 
voltage  in  (db)  and  the  third  axis  is  the  frequency  in  KHZ.  As 
the  frequency  is  changed  the  amplitude  will  rise  until  it  reaches 
a  maximum  and  then  decays  as  the  frequency  increases.  Figure  8  is 
a  right-hand  view  and  Figure  9  is  a  left-hand  view  of  the  same 
data. 

In  order  to  compare  the  theoretical  dispersion  curve 
versus  the  experimental  one,  the  data  for  the  experimental 
transducer  was  normalized  to  325  KHZ  so  that  the  comparison  could 
be  achieved.  The  normalized  dispersion  curve  is  plotted  for  the 
continuous  plate  in  Figure  10.  The  abscissa  is  the  steering  angle 
in  degrees  and  the  ordinate  is  the  frequency  in  KHZ.  The  points 
obtained  from  measurements  form  a  well-defined  curve  with  very 
little  scatter.  Although  the  experimental  points  are  similar  in 
shape  to  the  theory,  the  magnitude  of  dispersion  effect  of  the 
experimental  curve  is  significantly  less.  The  difference  between 
the  experimental  and  theoretical  curves  could  have  been  due  to 
limitations  in  the  model.  Water  loading  of  the  transducer. 


ignoring  boundary  impedance,  and  the  assumption  of  an  infinite 
plate  are  all  factors  that  could  influence  the  dispersion. 

3.2.2  Grating  Lobe 

In  Chapter  II  it  was  shown  theoretically  that  because  of 
the  spatial  filtering  in  the  continuous  plate  the  grating  lobe  in 
a  continuous  plate  was  reduced  considerably.  For  the  discrete 
segmented  plate  the  grating  lobe  occurs  when  the  spacing  between 
the  segments  is  greater  than  one-half  wavelength  when  the  beam  is 
steered  off  normal;  i.e.: 

KD(sin9  -  sin0o)  =  2n*  (n  =1,2,...) 

where  K  is  the  wavenumber,  0  is  segment  spacing,  9  is  the  physical 
angle  and  0Q  is  the  steering  angle. 


This  results  in: 


When  the  phase  between  the  segments  reaches  2r,  the 
pattern  repeats  itself,  and  this  happens  when  we  have  some  integer 
multiple  of  spatial  harmonic  of  the  wavenumber.  Both  the 
theoretical  and  experimental  data  for  the  discrete  segmented  plate 
show  such  an  effect.  The  theoretical  and  experimental  data  for 
10°  steering,  Z\  spacing,  four  bit  quantization  and  uniform 
shading  is  shown  in  Figures  11  and  12,  respectively.  The  first 
grating  lobe  in  both  cases  occurs  at  about  -27°,  as  predicted  by 
the  theory.  The  other  two  grating  lobes  can  also  be  seen  at  about 
+42°  and  -66°  in  Figures  11  and  12.  The  side-lobe  structure  of 
both  the  theoretical  and  experimental  pattern  seem  to  match  each 
other  rather  well. 

In  the  continuous  plate,  because  of  the  spatial  filtering, 
the  grating  lobe  was  observed  to  be  reduced  considerably,  as 
predicted  by  the  theory.  The  spatial  filtering  effect  in  the 
continuous  plate  was  measured  and  compared  against  the  theoretical 
data.  To  obtain  the  experimental  spatial  filtering  curve,  a 
particular  reference  steering  angle  was  selected  and  the  operating 
frequency  was  set  at  the  corresponding  resonant  frequency.  The 
amplitude  was  then  measured  for  different  electrical  drive 
wavenumbers.  The  normalized  experimental  curve  is  plotted  in 
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Theoretical  beam  pattern,  10°  steering,  2X  spacing,  4  bit 
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Figure  13  along  with  the  theoretical  spatial  filtering  curve.  The 
theoretical  spatial  curve  was  truncated  to  unity  to  achieve  a 
normal ization  equivalent  to  measurement.  This  plot  also  includes 
an  adjusted  wavenumber  normalization  which  was  a  best  match  to  the 
experimental  curve.  The  wavenumber  adjustment  necessary  to  match 
this  was  approximately  10%.  This  plot  Is  for  10°  steering.  The 
abscissa  is  the  normalized  wavenumber  and  the  ordinate  is  the 
normalized  amplitude. 

Figures  14  and  15  show  the  theoretical  beam  pattern  for  a 
10°  steering  angle,  2\  spacing,  and  uniform  shading,  with  no 
spatial  filtering  and  with  spatial  filtering,  respectively.  The 
first  grating  lobe  from  the  main  beam  has  been  reduced  by  16  db. 
The  other  two  grating  lobes  are  shown  to  have  an  improvement  of 
about  18  db.  The  experimental  beam  pattern  for  a  continuous  plate 
with  the  same  characteristics  is  shown  in  Figure  16.  The  first 
grating  lobe  is  down  at  approximately  -23  db,  which  agrees  well 
with  the  theory.  Although  the  other  two  grating  lobes  are  shown 
to  also  be  reduced,  the  reduction  is  about  15  db  versus  a 
theoretical  20  db. 
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Fig.  13  Theoretical  and  experimental  acoustic  amplitude  distribution 
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Theoretical  beam  pattern,  10°  steering,  2A,  using  FFT, 
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Experimental  beam  pattern,  10°  steering,  2x 
quantization 


3.2.3  Mirror  Lobe 

The  reflection  from  the  edge  boundary  of  the  continuous 
plate  that  was  discussed  in  Chapter  I  will  cause  a  standing  wave 
to  be  formed.  This  results  in  a  mirror  lobe  when  the  beam  is 
steered  off  normal.  An  active  approach  to  suppression  of  the 
mirror  lobe  was  used  and  gave  good  results.  Since  the  phase  for 
each  segment  could  be  set  independently,  it  was  possible  to  drive 
the  last  three  segments  in  the  plate  out  of  phase  with  the 
original  distribution.  This  acts  to  cancel  the  amplitude  of  the 
wave  by  the  time  it  reaches  the  edge  boundary  and  in  effect 
reduced  the  mirror  lobe  signi ficantly.  The  number  of  segments 
used  to  reduce  the  mirror  lobe  was  determined  experimentally. 

Figure  17  shows  the  12°  steered  beam  for  a  continuous 
plate  with  one-half  wavelength  spacing  between  thesegments.  The 
mirror  lobe  can  be  seen  to  have  a  high  level  nearly  comparable  to 
the  main  lobe  at  about  -12°.  Figure  18  shows  the  same  beam 
pattern  except  the  last  three  segments  were  driven  out  of  phase. 
The  mirror  lobe  can  be  seen  to  be  suppressed  considerably.  The 
suppression  is  about  12  db.  Figures  19  and  20  show  the  3-D  view 
of  the  fifteen0  beam  for  one-half  wavelength  spacing  for  the 
right-hand  view  and  the  left-hand  view,  respectively.  The  figures 
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show  how  the  mirror  lobe  moves  as  the  frequency  changes.  This  is 
an  indication  of  the  frequency  dependency  of  the  mirror  lobe  in 
the  continuous  plate. 

An  unexpected  low  level  (-15  db  or  less)  lobe  in  the 
measurements  of  the  beam  pattern  for  the  continuous  plate  was 
observed  at  approximately  -42°.  This  was  probably  generated 
because  of  the  shear  modes  in  the  plate  which  were  weakly  excited 
by  the  electrical  drive.  With  this  shear  mode  there  is  an 
associated  shear  velocity.  By  Snell's  law  the  angle  of  the 
radiated  beam  can  be  written  as: 

sine  * 

S 

where  Cw  is  the  speed  of  sound  in  the  water  and  Cs  is  the  shear 
velocity  in  the  plate.  In  the  continuous  piezoceramic  plate  this 
was  given  by  CHANNEL  INDUSTRIES  to  be  2250  meters/sec.  Although 
the  location  of  the  shear  lobe  was  independent  of  both  steering 
angle  and  frequency,  the  amplitude  of  the  radiated  lobe  had 
dependency  on  both  frequency  and  steering  angle. 

Figure  19  shows  a  right-hand  view  and  Figure  20  shows  a 
left-hand  view  for  a  15°  steered  beam  and  one-half  wavelength 


spacing.  The  abscissa  is  the  angle  in  degrees,  the  ordinate  is 
the  normalized  voltage  in  db  and  the  third  axis  Is  the  frequency 
in  KHZ.  The  effect  of  frequency  can  be  seen  to  change  the 
amplitude  of  the  shear  lobe  at  -42°s  in  the  beam  pattern  of  the 
continuous  plate.  The  effect  of  steering  angle  on  the  shear  lobe 
is  illustrated  by  Figure  21,  which  shows  a  right-hand  view,  and 
Figure  22,  which  shows  a  left-hand  view  for  a  0°  beam.  The 
amplitude  of  the  lobe  is  about  -24  db  for  0°  steering  versus  the 
-16  db  for  15°  steering. 


3.2.4  Quantization  Effect  on  the  Lobes 


of  the  Discrete  and  Continuous  Plate 


The  effect  of  four,  three  and  two  bit  quantized  driving 
signals  was  investigated  for  the  discrete  and  continuous  plates. 

Figure  23a,  23b  and  23c  shows  a  10°  steered  beam  for  a 
discrete  plate  at  1/2X  spacing,  uniform  shading  for  four,  three 
and  two  bit  quantization,  respectively.  Because  of  the  finer 
phases  of  the  4  bit  quantization,  the  side-lobe  structure  for  this 
quantization  was  best  in  comparison  with  the  other  two 
quantizations.  This  effect  can  be  seen  very  well  at  -23°,  where 
the  side-lobe  is  about  -20,  -15,  -14  db  for  the  four,  three  and 
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Fig.  23c  Beam  pattern,  \  2  bit  quantization,  discrete  plate 
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two  bit,  respectively.  The  first  side-lobe  of  the  main  beam  for 
the  four  bit  is  -9.5  db  versus  the  side-lobe  of  about  -8.5  db  and 
-8  db  for  the  three  and  two  bit,  respectively.  The  better 
performance  of  the  four  bit  quantization  can  also  be  seen  at  -50° 
where  the  side-lobe  is  at  -24  db  versus  the  side-lobes  of  -20  db 
and  -16  db  for  the  three  and  two  bit  cases,  respectively.  The 
effect  of  quantization  ranges  (2-6  db)  at  -23°  and  (1-2  db)  for 
the  first  side-lobe. 

Figure  24a,  24b  and  24c  shows  a  5°  steered  beam 
for  1/2X  spacing  in  a  continuous  plate  for  four,  three  and  two 
bit,  respectively.  All  side-lobes  for  the  continuous  plate  are 
lower  than  those  of  the  discrete  plate  and  the  effect  of  phase 
quantization  is  also  smaller  for  the  continuous  plate.  The  first 
side-lobe  is  about  -15  db  for  all  the  quantizations.  The  four  bit 
has  a  side-lobe  of  about  -19  db  at  a  28°  angle.  The  three  and  two 
bit  shows  the  same  side-lobe  at  about  -14  db  and  -17  db, 
respectively,  which  is  not  too  much  improvement  at  that  particular 
angle.  But  at  -12°  the  side-lobe  for  the  two  bit  is  about  -16  db 
versus  -24  db  of  the  three  bit  quantization.  This  is  rather  a 
good  improvement  in  the  side-lobe  structure  of  a  continuous  plate. 


FIVE  DEC-  ETEE^I 


,  4  bit  quantization,  continuous  plate 
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The  grating  lobe  In  the  discrete  plate  also  shows  a 
dependency  on  the  quantization.  Figure  25a,  25b  and  25c  shows  a 
5°  steered  beam  for  2X  spacing  for  four,  three  and  two  bit, 
respectively.  The  first  grating  lobe  for  four  bit  quantization  is 
about  -14  db.  The  three  and  two  bit  also  show  that  they  have 
grating  lobes  of  about  -12  db  and  -10  db,  respectively,  at  about 
-28°  angle.  The  four  bit  quantization  is  shown  to  have  a 
considerably  reduced  grating  lobe  compared  to  the  others. 

For  a  continuous  plate,  the  grating  lobe  was  considerably 
reduced  because  of  the  spatial  filtering  effect.  The  effect  of 
coarser  quantization  is  also  shown  to  be  less  than  for  the 
discrete  plate.  Figure  26a,  26b  and  26c  shows  a  5°  steered  beam 
for  2X  spacing.  The  grating  lobe  is  about  -30,  -24,  -24  for  four, 
three  and  two  bit,  respectively. 

For  the  continuous  plate  the  mirror  lobe  was  suppressed 
more  effectively  with  four  bit  quantization  because  of  the  finer 
phase  adjustment  that  was  available  to  drive  the  last  three 
segments  for  cancellation.  Figure  27a,  27b  and  27c  shows  a  10° 
steering  for  IX  spacing,  uniform  shading  for  four,  three  and  two 
bit  quantization,  respectively.  The  four  bit  quantization  effect 
on  the  mirror  lobe  is  about  -22  db  versus  the  -20  and  -19.5  db  for 
the  three  and  two  bit,  respecti vely. 
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4.0  Conclusion 

The  purpose  of  this  thesis  was  the  study  of  the  radiating 
properties  of  a  segmented  continuous  piezoelectric  plate 
transducer  in  comparison  with  a  discrete  segmented  plate. 
Continuous  plates  have  been  shown  to  have  both  advantages  and 
disadvantages.  The  major  advantages  of  the  continuous  plates  are 
the  spatial  filtering  effect  and  ease  of  fabrication.  The 
disadvantages  are  spurious  shear  lobes  and  mirror  lobes,  and  a 
narrow  operating  frequency  bandwidth. 

The  spatial  filtering  effect  of  the  continuous  plate  had  a 
very  noticeable  impact  on  the  grating  lobes.  This  effect  was 
evaluated  both  theoretically  and  experimentally.  The  theoretical 
analysis  was  carried  out  by  computing  the  response  of  the  plate  in 
wavenumber  space  and  using  that  response  as  a  weighting  function 
on  the  wavenumber  spectrum  of  the  impressed  driving  voltage. 
Experimental  measurements  confirmed  the  theoretical  predictions. 
The  first  grating  lobe  for  the  continuous  plate  was  down  to  -29  db 
compared  to  -4  db  for  the  discrete  plate. 

The  spatial  filtering  effect  of  the  continuous 
piezoelectric  plate  transducer  served  to  smooth  the  side-lobe 


structure  significantly.  The  side-lobe  performance  of  the 
continuous  plate  was  lower  than  that  of  the  discrete  plate.  The 
first  side  lobe  for  the  continuous  plate  was  down  -12  db  compared 
to  -8  db  for  the  discrete  plate. 

The  spatial  spectrum  of  quantized  drive  for  different 
phasing  was  also  investigated  for  four,  three  and  two  bit.  The 
four  bit  quantization  gave  the  best  performance  compared  to  the 
two  bit  quantization  for  both  continuous  and  discrete  plate.  The 
first  side  lobe  for  the  four  bit  quantization  in  a  continuous 
plate  was  down  -12  db  compared  to  -10  db  for  the  two  bit 
quantization  for  2X  spacing  and  10°  steering,  a  degradation  of 
2  db  due  to  quantization.  The  first  side  lobe  for  the  four  bit 
quantization  of  the  discrete  plate  was  down  -9  db  compared  to 
-3  db  of  the  two  bit  quantization  for  2\  spacing  and  10°  steering, 
showing  a  greater  degradation  of  6  db. 

In  the  continuous  plate  the  mirror  and  shear  lobes  were 
shown  in  Chapter  II  to  be  of  major  concern.  The  mirror  lobe  was 
generated  by  the  reflection  of  the  travelling  wave  from  the  edge 
boundary  which  created  a  standing  wave.  The  mirror  lobe  was 
experimentally  cancelled  by  actively  driving  the  last  three 
segments  out  of  phase,  which  reduced  the  mirror  lobe  by  about 


18  db.  Shear  lobes  were  experimentally  observed  to  be  at 
approximately  42°.  These  were  symmetrical  shear  waves  which  were 
weakly  coupled  to  the  electrical  drive. 

A  suggestion  which  could  be  made  here  is  to  passively 
suppress  both  the  mirror  lobe  and  shear  lobe  by  using  a  proper 
damping  at  the  edges  of  the  continuous  plate  to  eliminate  a 
standing  shear  wave  or  normal  mode  which  might  be  generated  in  the 
continuous  plate. 

A  very  important  physical  advantage  of  the  continuous 
plate  over  the  discrete  plate  may  be  in  the  ease  of  fabrication, 
especially  for  high  frequencies.  At  high  frequencies, 
wavelengths  are  small  and  fabrication  of  a  multi -element  array  is 
difficult,  particularly  for  imaging  sonars  working  in  the 
megahertz  region  with  apertures  on  the  order  of  100  wavelengths. 
For  narrow  scanning  apertures  the  conducting  electrode  of  the 
continuous  plate  can  be  cut  at  rather  large  spacing  (greater  than 
1/2  wavelength)  without  worrying  about  the  grating  lobes,  since 
they  will  be  reduced  by  the  spatial  filtering  effect  of  the 
continuous  plate.  The  necessity  of  dicing  the  fragile  ceramic 
plate  may  also  be  circumvented  by  merely  using  individual 
electrodes  created  by  silk  screening  the  electrode,  pattern 
initially  or  created  by  etching  a  continuous  electrode  surface. 
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